function P = BB98d_predict(T_a,dT_a_dt,x)

    clear('P')
    P.a  = x(1);
    P.b  = x(2);
    P.c  = x(3);
    if numel(x) > 3
        P.amp_a  = x(4);
        P.phi_a  = x(5);
    end
    if numel(x) > 5
        P.amp_b  = x(6);
        P.phi_b  = x(7);
        P.amp_c  = x(8);
        P.phi_c  = x(9);
    end
    P    = BB98d_generate_data(T_a,dT_a_dt,P);
    
end
